A resolution to doping asymmetry puzzle in high T c cuprates 
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We present a microscopic model for 'electron doped' Nd2-a:Ce :I :Cu04 family and offer a resolution 
to a long standing doping asymmetry puzzle. Here, i) Ce atoms do not dope free electrons, instead 
a Ce atom effectively quenches a Cu 2+ spin moment at an adjacent site, at an energy scale > 
superexchange J of CuC>2 plane and 'site dilutes' the Mott insulator, ii) effective chemical pressure, 
caused by increased Ce substitution, induces a first order Mott insulator to superconductor transition 
in the Cu02 plane. We predict, i) equal number of +e and -e carriers in metallic state and ii) a line 
of first order transition ending at a critical point in normal state. Phenomenology gets organized. 
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Following the seminal discovery of superconductivity in 
La2- x Ba x Cu02 |l|, the idea of resonating valence bond 
mechanism of superconductivity in 'hole doped Mott in- 
sulator' was bornQ. In this theory 0, based on a single 
band model, electron doping of Mott insulator should 
equally well lead to superconductivity and one expects 
some kind of symmetry between electron and hole dop- 
ing. When a so called 'electron doped' superconductor 
Nd2- a; Ce a ;Cu04 was discovered^ there was a surprise. 
At the level of phase diagram (figure la) there was a strik- 
ing asymmetry and significant differences. A large body 
of experiments [4-23] devoted to a study of this asymme- 
try exists now. In spite of several theoretical |24| at- 
tempts, our understanding of this asymmetry is not sat- 
isfactory. 

Antiferromagnetism in Nd2- x Ce ;E Cu04 family extends 
to a doping of nearly 12 %, in contrast to a 2% doping 
in LSCO (figure 1). Superconductivity appears abruptly 
at a first order transition with a maximum T c which de- 
creases on further Ce doping^, 0, 0|. Hall and other 
transport measurements show anomalies 0, B- different 
from that shown by hole doped cases. 

The large asymmetry naturally points to some missing 
features in the current modeling of superconductivity in 
Nd2- a; Ce a ;Cu04 family. In the present article we identify 
these features and propose a new model, where strong 
electron correlations continue to play major role. The 
missing features follow from some known experimental 
results and some physical arguments. Our model is con- 
sistent with the large body of existing phenomenology. 
It also makes some predictions. 

We argue that root cause of doping asymmetry lies 
in the T'-structure of Nd2^ 2; Ce a ;Cu04 family and am- 
bivalent (sic) character of Ce valence. T" differs from 
T, the structure of La2Cu04 : apical oxygen of Cu04 
octahedra (T-structure) are displaced in the ab-plane 
and fall in line with the oxygens of CuC>2 layers along 
the c-axis. A remarkable consequence of missing apical 
oxygen in T'-structure seems to be a nearly 1 eV re- 
duction in Mott-Hubbard gap in Nd2- a; Ce a: Cu04 (com- 
pared to La2Cu04 ), as seen in recent x-ray scattering 
experiment [13I ■ This makes Nd2Cu04 a fragile Mott in- 
sulator, close to a Mott insulator to superconductor tran- 
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FIG. 1: a) Experimental phase diagram of electron and 
hole doped cuprates, reproduced from refff. b) dop- 
ing dependence of Neel temperature for electron doped 
Pr2-iCexCu04 , oxygen doped La2Cu04+i and Zn doped 
La2Cui_ y Zn y 02, reproduced from reference|(J. 



sition point. The unstable tetra valent state of Ce is 
likely to play a key role in quenching a Cu 2+ spin mo- 
ment, and also generate required chemical pressure for 
a Mott insulator to superconductor transition, that we 
propose below. 

In what follows, we define our microscopic model and 
discuss its solution. We also indicate how existing quan- 
titative theoretical understanding on the hole doped case 
can be transfered to the Nd2-a;Cea;Cu04 family, because 
of some special feature of our model. 

Hypothesis leading to our model are: A) CuC-2 planes 
are not electron doped by Ce reservoir, instead every Ce 
effectively quenches a Cu 2+ spin moment of an adjacent 
site, over an energy scale > superexchange J of CuC>2 
plane and 'site dilutes' the Mott insulator. B) chemical 
pressure arising from increased Ce substitution, induces 
a first order Mott insulator to superconductor transition 
in the CUO2 plane, resulting in a 'self doping' of Mott 
insulator with equal number of +e and -e carriers. 

Our model is 'a site diluted Heisenbcrg model' for the 
insulating phase, which exists in the experimental up to 
a value of x c ~ 0.12. For x > x c , we have a conductor, 
described by a 'site diluted 2-species t-J model', which 
we elaborat later. In both cases, number of 'missing Cu 
sites' or 'diluted sites' is Nx, where N is the number of 
sites and x is the Ce density. 
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It should be mentioned that effective 'site dilution' as 
an experimental consequence of Ce doping (in the insu- 
lating antiferromagnetic order regime) is known in the 
liter ature[(|, including a recent theoretical fit of their 
neutron scattering result to the site diluted Heisenberg 
model by Mang et al.0. Results, which is behind the site 
dilution hypothesis are: i) the Neel temperature Tn has a 
slope which is the same as the slope for a well established 
site dilution situation namely La2Cui_ x Zn x 02 , and ii) 
the long range antiferromagnetic order continues to be a 
commensurate (tt, 7t) order. These two results are qual- 
itatively different from the hole doped case, where Tn 
is suppressed completely by about 2 % doping and any 
magnetic quasi long range order is incommensurate and 
not at (7r,7r). We take site dilution as a phenomenologi- 
cal input to our modeling and do not discuss mechanism 
and theory behind it. As far the theory of site diluted 
2D Heisenberg model, we refer to existing work in the 
literature. 

Two sceneries are possible in metallic state (x > x c ) of 
Nd2- a; Ce a ;Cu04 : i) Ce does not quench a Cu 2+ moment 
(i.e., no site dilution), instead becomes Ce 4+ and elec- 
tron dopes the CuC-2 layer, and ii) spin quenching (site 
dilution) continues and we reach metallic state through 
a Mott transition, induced by chemical pressure. Recent 
observation of Kondo like effect [l2^ in Pr2- x Ce x Cu04 , 
according to us, indicates that Ce continues to be in- 
volved in spin quenching in the metallic state. Further 
there is a strong indication of first order insulator to 
metal transition 13J, reminiscent of of ET-salts (organic 
superconductor), where pressure drives a Mott insulator 
to superconductor phase transition^. Thus we hypoth- 
esize that, based on the above and overall phenomeno- 
logical grounds, our system undergoes a transition from 
a (site diluted) Mott insulator to (site diluted) supercon- 
ductor. The pressure in the present case is an effective 
chemical pressure, arising from reorganization of electron 
wave functions and hopping parameters, arising from a 
strong electronic coupling to the cerium subsystem. 

What do we know about pressure driven Mott insula- 
tor to superconductor transition ? The present author- 
recently predicted 26] such a possibility in cuprates and 
other systems; Zhang |27j has an independent and some 
what different approach. In our mechanism long range 
coulomb interaction, originally invoked by Mott plays a 
key role in establishing a strong first order transition. 
Our key new input is to view the metallic state close to 
the Mott transition as a 'projected metal', a 'self doped' 
Mott insulator, where superexchange survives in the con- 
ducting state. Within a simple tight binding model, the 
transition is viewed as a spontaneous creation of a small 
but finite and equal density of -e (doubly occupied orbital 
or 'doublon') and +e (empty orbital or 'holon') charge 
carriers, whose individual number is approximately con- 
served at low energies. What is important is that the 
superexchange survives in this conducting state - this is 



a memory of the Mott insulator. The concentration of 
self doped carriers 2y (= y + y for e~ and e + carriers) is 
decided by long range coulomb interaction, local dielec- 
tric constant and band parameters. A projected metal is 
meaningful as long as 2y << 0.5 (fig. 3b). Here 0.5 is den- 
sity of doubly occupied and empty orbitals in one orbital 
tight binding model of free electrons at half filling. 

In our work|26j, proposing a mechanism for Mott in- 
sulator to superconductor transition under pressure, we 
introduced a 2-species t-J model, to take care of the 
projected character of the self-doped state, and dis- 
cussed superconductivity. While the insulating phase of 
Nd2-a;Ce x Cu04 is described by a Heisenberg model on a 
diluted lattice, the conducting state is described by our 
2-species t-J model on a diluted lattice: 

H 2 u ~ - XI %4r s ja( e ; e ] - d i d ]) + h - c - 

ij 

including local constraint d\d i + e[e i + J2a s la s ia = 1- 
Summation over sites excludes 'diluted' sites. Here we 
have used the slave boson representation for an electron 
c ia = e]s icr +(7d i sl s , with e's, d's and s's representing zero 
occupancy (+e), double occupancy (-e) and single occu- 
pancy (neutral). The first term of equation 1, represents 
holon hopping and doublon hopping, the second term 
the superexchange among neutral sites. The conserved 
number of singly occupied sites is fixed through a global 
constraint, Ylia- s \<y s ia = (1 — 2y)A; this in turn also fixes 
the holon and doublon numbers to be yN. As mentioned 
earlier, y is a self consistent parameter that depends on 
the chemical pressure, in our case on Ce concentration 
x. The physics of long range coulomb interaction that 
caused a first order Mott transition is contained in our 2- 
species t-J model, parametrically through the self doping 
density y. 

How superconductivity arises in 2-species t-J model ? 
We have shown earlier |26| that for the case of bipartite 
lattice a 2-species t-J model can be exactly mapped onto 
the regular (single species) t-J model. This mapping fol- 
lows once we recognize that d's and e's are hard core 
bosons and they also commute at different sites. The 
transformation di — > implements the exact mapping. 
Here e,± on the A and B sublattices of the square lattice. 
As site dilution preserves the bipartite character of the 
square lattice, the mapping survives and we get a regular 
t-J model on the same diluted lattice. 

Presence of (a small) next nearest neighbor hopping f 
makes the transformation inexact. However, in such a 
case one can study the 2 species t-J model directly using 
various known approximate methods. 

What is the effect of site dilution on the symmetry 
of superconducting order parameter. Site dilution is a 
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Predicted line of first orderr 
metal-insulator (Mott) transition 
ending at a critical point 



Mott Insulator 1 
(site diluted) 



AFM Order 
commensurate (71,71 ) 



Self doped 
Mott Insulator 
(site diluted) 




Ce-doping x (Chemical pressure) — 

FIG. 2: Our proposed phase diagram. A key prediction is 
the presence of a line of first order metal-insulator (Mott) 
transition ending at a critical point 



strong perturbation. One would expect, based on An- 
derson's theorem on dirty superconductors, that T c of d- 
wave superconductivity will be very sensitive to the site 
dilution disorder. Here, something remarkable happens 
and makes site dilution disorder less effective in suppress- 
ing d-wave superconductivity. We will not go into the 
details, except to point out that this happens within our 
RVB mean field analysis for the following reason. A gen- 
eral scatterer will cause phase randomization as cooper 
pairs scatter in k-space, because of the change in sign 
(arising from d x 2_ y 2 symmetry). However, our scatterer 
is a site dilution and phase randomization does not occur, 
as two sublattice character is preserved by site dilution. 
That is, in real space, the d-wave order parameter fits 
into a site diluted 2D square lattice without any frustra- 
tion. 

The above is consistent with what is experimentally 
known about symmetry of the superconducting order 
parameter 0| and also the scale of superconducting 
T c for the Nd2_ :c Ce a; Cu04 family. It should be pointed 
out that a d x i_ y i to extended-S transition suspected^] 
in some experiments is allowed by our theory, as site 
dilution will have different quantitative effects on an 
extended-S solution and d x 2_ y 2 solution and might stab- 
lise s-wave at higher dopings. 

In the light of our model the following experimental 
results get organized and are less puzzling, and also some 
predictions follow. We summarize them below: 

i) In a recent work, Tsukada et al..|lfij found super- 
conductivity in La2CuC>4 in T'-structure, for trivalent 
rare earth substitution, Laj^R^CuC^ ; R — Lu and 
Y, among others. This is an unexpected and remark- 
able finding. A well known stable trivalent state of Lu 
and Y makes doping of CuC>2 plane unlikely. We suggest 
that La2CuC>4 in T'-structure has an inherent enhanced 
chemical pressure, which is further increased by Lu or 
Y substitution leading to Mott insulator superconductor 
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FIG. 3: a) Mott-Hubbard or Charge gap E c as a function of 
chemical pressure. Schematically, three experimentally stud- 
ied systems have different charge gaps as shown, b) Self-doped 
carrier density y as a function of chemical pressure (related to 
Ce doping x). Self doping looses its meaning as y approaches 
0.25, and we leave a projected metal and enter a fermi liquid 
like state 



transition without external doping, via our mechanism. 

A recent inelastic X-ray scattering study 0] gives 
a Mott Hubbard gap of Nd2- a; Ce a ;Cu04 smaller than 
LSCO by about 1 eV. This is an experimental evidence 
for different quantum chemistry in T and T' structure 
and that Nd2- a; Ce a ;Cu04 is under higher chemical pres- 
sure and close to a Mott transition point compared to 
La2CuC>4 (figure 3a). It will be important to measure the 
Mott Hubbard gap in La2CuC>4 in T'-structure, which 
according our suggestion will also be close to the Mott 
transition point. 

ii) The shape of the superconductor-antiferromagnetic 
phase boundary and various experiments 13] indicate a 
strong first order character of the transition and phase 
separation. Our hypothesis B not only explains it in 
a most natural fashion but also makes the following 
important prediction. On very general thermodynamic 
grounds, the first order transition (as it is driven by 
coulomb interaction among charges) will continue beyond 
the antiferromagnetic phase boundary into the paramag- 
netic normal state and end at a critical point, as shown 
in figure 2. This is experimentally established in organic 
superconductors, where pressure drives a Mott insulator 
to superconductor transition. 

When the critical concentration x c is small, disorder 
will have less effect on the first order line and the end 
critical point. As smaller critical x c is realized in an- 
nealed Pr2_ a; Ce a ;Cu04 and La2Cu04 systems, it may be 
easier to confirm our prediction of first order transition 
line and critical point in the normal state. Nuclear and 
muon spin resonance will be good local probes to study 
the above. Further, as the first order phase transition is 
charge driven, there will be a volume change across the 
phase boundary, due to coupling of the charges with the 
lattice. This can be detected by sound velocity, elastic 
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FIG. 4: STM spectra, reproduced from reference |2(j|). for 
various tip to surface distance X. According to our theory, 
Nd2- a; Ce :E Cu04 , which is self-doped with equal density of +e 
and -e carriers, has a symmetric spectra. Hole doped YBCO, 
which has only +e carriers has an asymmetric spectra, x-axis 
scale is meV 



constant studies. 

iii) In a recent paper, Anderson and Ong|l°J] have sug- 
gested that tunneling conductance should be intrinsically 
asymmetric, in either electron or hole doped Mott insu- 
lators. We have applied their analysis to our 2-species 
t-J model and find a symmetric tunneling spectra, as ex- 
pected on symmetry grounds. There is some asymmetry 
to the extent of particle-hole asymmetry (arising from 
t' etc.). Recent results show |2fl 12 1| a high symmetry in 
the tunnelling conductance. In figure 4, the high symme- 
try in N do-y Ce^CuOa compared to the case of YBCO of 
reference [20j is exhibited; it is consistent with our pro- 
posal. To substantiate our model further, it will be in- 
teresting to perform more study on symmetry aspect of 
tunneling in Nd2- a; Ce a ;Cu04 , as well as other systems 
with T' structure. 

iv) Earlier electrical transport studies have invoked^ 
the presence of two species of charges to explain certain 
anomalous results. It will be interesting to look for our 
prediction of equal density of hole and electron carriers in 
the conducting state of Nd2-^Ce :c Cu04 family and other 
T' structures, using other experiments. 

v) There has been extensive ARPESj^ results on 
fermi surface measurements in Nd2- £! ;Ce a ;Cu04 family. 
ARPES finds a Luttinger volume of ( 1 ~^) JV k-points. 
This is consistent with our model - our conducting half 
filled band, in view of site dilution, has same number of 
effective k-points as a full lattice with an electron doping 
density of x.) 

vi) Recent NMR results and other experiments ^23[ in- 
dicate a 'over doped' fermi liquid like behavior in the 
normal state, We believe that this is due to the fact that 
our effective carrier density 2y is beyond the standard 
optimal doping density of 0.15, making it effectively over 
doped. As we discussed earlier, Ce density x is different 
from self doping density 2y; and 2y is determined by band 
parameters, dielectric constant and coulomb interaction 
physics. 
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